Systematic study of Optical Feshbach Resonances in an ideal gas 
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Using a narrow intercombination line in alkaline earth atoms to mitigate large inelastic losses, we 
explore the Optical Feshbach Resonance (OFR) effect in an ultracold gas of bosonic 88 Sr. A system- 
atic measurement of three resonances allows precise determinations of the OFR strength and scaling 
law, in agreement with coupled-channels theory. Resonant enhancement of the complex scattering 
length leads to thermalization mediated by elastic and inelastic collisions in an otherwise ideal gas. 
OFR could be used to control atomic interactions with high spatial and temporal resolution. 
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The ability to control the strength of atomic interac- 
tions has led to explosive progress in the field of quan- 
tum gases for studies of few- and many-body quantum 
systems. This capability is brought about by magnetic 
field-induced Feshbach scattering resonances (MFR) Q, 
where both the magnitude and sign of low-energy atomic 
interactions can be varied by coupling free particles to 
a molecular state. MFR in ultracold alkali atoms have 
been used to realize novel few-body quantum states and 
study strongly correlated many-body systems and phase 
transitions |2[. However, magnetic tuning has limited 
current experiments to relatively slow time scales and low 
spatial resolution. Higher resolution could be achieved by 
controlling MFR optically 0]. 

Scattering resonances can also arise under the influence 
of laser light tuned near a photoassociation (PA) reso- 
nance 0] where free atom pairs are coupled to an excited 
molecular state 0, Q . This Optical Feshbach Resonance 
(OFR) is expected to enable new and powerful control 
with high spatial and temporal resolution. OFR has been 
studied in thermal 0] and degenerate [H 0| gases of Rb, 
but it was not found useful due to large photoassociative 
losses. Much narrower optical intercombination lines are 
available in alkaline earth atoms and are predicted to 
overcome this loss problem fioj ]. Independently, ultra- 
cold alkaline earth atoms have recently emerged to play 
leading roles for quantum metrology [LU - llla ] where preci- 
sion measurement and many-body quantum systems are 
combined to study new quantum phenomena [13, flij ]. 
Degenerate gases of alkaline earth atoms have recently 
become available [l6[ . Due to the lack of magnetic struc- 
ture in the ground state of these atoms, the OFR effect 
could become an important tool for controlling their in- 
teractions. OFR work on Yb [13, EH has been limited 
to studying the induced change in scattering phase shifts 
and PA rates. Dominant PA losses are evident in all of 
the OFR experiments listed above. Light-induced clastic 
collisions for thermalization were not observed. 

In this Letter, we study the OFR effect across mul- 
tiple resonances in a metastable molecular potential of 
88 Sr. The aim of this work is to test the practical appli- 



cability of OFR for engineering atomic interactions in the 
presence of loss, similar to the successful application of a 
decaying MFR For 88 Sr, OFR is predicted [ijj] to al- 
low changes in the scattering length by more than a factor 
of 100 with low losses by using large detunings (0(1O 5 ) 
linewidths) from the least-bound vibrational level poj ]. 
We tested this proposal and find experimentally that the 
existing isolated resonance model [6J only describes the 
experiment in the small detuning regime. Large detun- 
ings from a molecular resonance require a full coupled- 
channels description of the molecular response. Sup- 
ported by this new theory framework, we present a sys- 
tematic experimental study of the OFR-enhanced com- 
plex scattering lengths and demonstrate OFR-induced 
thermalization in an ultracold gas. 

Bosonic 88 Sr has an s-wave background scattering 
length <2b g = — 1.4(6)arj (HI, where ao is the Bohr ra- 
dius. The small |ob g | makes the sample effectively non- 
interacting and provides an ideal testing environment for 
OFR. FigureQJt shows the ground ( 1 So- 1 So g ) and low- 
est excited state ( 1 So- 3 Pi U ) molecular potentials of Sr 2 , 
which are coupled by a PA laser near the atomic tran- 
sition at A Q = 689 nm. The vibrational levels investi- 
gated are labelled by their quantum number n, counted 
as negative integers from the free particle threshold. For 
a given PA laser detuning from threshold, the Franck- 
Condon principle localizes the atom-light interaction in 
the vicinity of the Condon point Q . 

When detuning the PA laser across a vibrational res- 
onance, the s-wave scattering length shows a dispersive 
behavior, just as for a MFR. However, the finite lifetime 
of the excited molecular state leads to loss intrinsic to 
OFR. This process can be described [22[ akin to decay- 
ing MFR 0, [l9[ with a complex s-wave scattering length 
a(k) = a(k) — ib{k) that depends on the relative momen- 



tum hk and a PA line strength factor 
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called the optical length [lfj, |23| . Here, c is the speed of 
light, and £ op t scales linearly with PA intensity / and 
free-bound Franck-Condon factor |(n|.E)| 2 per unit colli- 
sion energy E = h 2 k 2 /(2fj,) at reduced mass [i = msr/2. 
In the isolated resonance approximation 0] the inelastic 
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FIG. 1: (a) Ground (blue) and excited (red) molecular po- 
tentials of Sr2. The inset shows the difference potential after 
subtracting the optical frequency. Horizontal lines indicate 
bound molecular states n in the excited potential. The free 
particle (bound state) radial wave function is indicated in blue 
(red), (b) Loss spectra for U n=-2, -3, and -4 for exposure 
time tpa = 200 ms and comparable mean density. / is scaled 
to keep F sc sufficiently small. The similarity of the spectra 
demonstrates the universal scaling with £ op t oc |(n|i?)| I. 



collision rate is [22 1 



K in (k) = 



lopt ~jn 



M (A + E/h) 2 /-/ 2 + [l + 2k^f^} 2 /4 



(1) 



where A is the laser detuning from molecular reso- 
nance 22j. We have accounted for extra molecular 
losses with 7 > j m = 2~f a , where j m is the linewidth 
of the molecular transition and 7 a = 2tt x 7.5 kHz 
is the atomic linewidth. Neglecting dbg for 88 Sr gives 
K c \(k) ~ 2k e ° pt J m Ki n (k) . The elastic-to-inelastic colli- 
sion ratio K e \/K- m becomes less favorable for smaller k. 

We load ^5xl0 4 atoms from a magneto-optical trap 
operating on the 1 So- 3 Pi intercombination transition 
into a crossed optical dipole trap formed by tilted hori- 
zontal (H) and vertical (V) beams (1064 nm), with trap 
depths ~15 fiK and ~7 /iK, respectively. The trapped 
sample shows a clear kinetic energy inhomogeneity be- 
tween the H and V axes (2-2.5 /xK vs. 3-4 /xK), due 
to the negligible ab g , consistent with a thermal distribu- 
tion energy-filtered by the trap potential. Typical in-trap 
cloud diameters are 45-55 /im. The PA beam intersects 
the trap with a waist of 41 /mi [22l |. 

A representative survey of PA resonances in the 1 So- 
3 Pi U potential is shown in Fig. [T|d. The PA laser with 
intensity /, adjusted to achieve similar £ opt for all spectra 
shown, interacts with the sample for tpa- Photon- atom 



scattering at rate T sc and subsequent radiation trapping 
set the maximum usable I for a given detuning from the 
atomic line [2^. In addition to the vibrational levels in- 
dicated in Fig. \T^l, the n=-l vibrational state exists at 
-0.4 MHz detuning from the threshold, which leads to a 
PA resonance with a very large line strength £ opt / 1 ■ 
The isolated resonance theory indicates that operating 
with a large I at 0(lO 5 7 a ) detuning from the n=-l state 
should allow modifications to a(k) of O(100ao) [ioj . This 
prediction relied on extrapolating the large line strength 
of the n=-l state across multiple intermediate PA reso- 
nances. However, with / up to 1 kW/cm 2 and detunings 
up to -1.5 GHz, we did not observe any effects due to 
elastic collisions. 

The discrepancy between theory and experiment stim- 
ulated a coupled-channels treatment of an atomic colli- 
sion in a radiation field that properly switches between 
the short range molecular states and two field-dressed 
separated atoms 22, 24, [25|. In the coupled-channels 
theory, the two coupled excited potentials (0 U , l u ) have 
the form of Ref. [20j, with an added imaginary term 
—ih'y m /2. The ground state potential uses the dispersion 
coefficients of Ref. [26| , has a scattering length of — 1.4 ao, 
and reproduces the bound state data of Ref. [2l| to better 
than 0.4%. Coupled-channels calculations do not assume 
isolated resonances, and all 0„ and l u molecular eigen- 
states emerge from the calculation as interfering, decay- 
ing scattering resonances [lj]. 

FiguresEt and d show that the coupled-channels model 
reproduces the isolated resonance expressions [l], 0] for 
a(k) and the rate constants as long as A is small com- 
pared to the spacing between molecular levels. How- 
ever, the coupled-channels K c \ returns to its background 
value K e i in between resonances regardless of their rel- 
ative strengths. The dotted line indicates (abg) a t 
E/k B = 4 fiK in Fig. Hi (Fig. Hi). These calculations 
show that each molecular line behaves as an isolated res- 
onance near its line center. For detunings comparable 
to the molecular level spacing, the isolated resonance ex- 
pressions cannot be used. 

At intermediate detunings, |A| 3> 7(1 + 2k£ opt ^), 
a(k) can be written in the standard form for MFR [22|, 



lim a{k) = ab g 

k~ >0 
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A 



2 A 2 



(2) 



To obtain a meaningful 
-opt 7m/ A needs to be suf- 



where w = — 4pt7m/abg- 
change in scattering length, 
ficiently large, and the imaginary part b = ^£ opt j m J / A' 2 
needs to be sufficiently small. Since K ln ~ (2 x 
10~ 12 cm 3 /s) (b/ao), for a density of p = 10 12 cm~ 3 and 
b = O.lao, K m p = r sc for I = 53 W/cm 2 assumed for 
Figs. [2j:,d. Thus, the calculations predict that changes 
in the scattering length of order 10 ao ^ \ a hg\ should 
be possible with O(100 j m ) detunings on timescales of 
200 ms. 
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FIG. 2: (a) Coupled-channels calculations of K e \ and .Kin at 
E/ks = 4 /iK and I = 10 W/cm 2 , versus PA laser detuning 
from atomic resonance. Each resonance peak is labeled by 
its electronic symmetry U or l u and n. Between resonances, 
K- m is only approximate, (b) Ratio of thermally-averaged rate 
constants at 2 /iK for A/7 m near U (— 2). / = 44 mW/cm 2 
gives the same ^ op t = 360ao as for the conditions in Fig. [4)3. 
(c) Zero energy limit of a(k) and b(k) for detuning near the 
U (— 4) feature, (d) Same calculation as in (c) at large detun- 
ing. Here, the isolated resonance results (solid lines) agree 
with the coupled-channels theory (circles). 



To investigate the utility of OFR, we systematically 
characterized three different resonances and determined 
their universal scaling. Because K c \/ K m oc k£ opt , inelas- 
tic collisions dominate the dynamics of the sample for 
small £ opt < (2(/c))- 1 



= | y/n/ T), where the an- 
gled brackets indicate a fc-average at temperature T, and 
ks is the Boltzmann constant. In this regime, the result 
of scanning the PA laser across resonance is a loss feature 
that shows no dependence on elastic collision processes. 

A typical PA loss feature for small £ opt is shown in 
Fig. [3^, where the final atom number after applica- 
tion of PA light is shown with respect to PA detuning 
from 1 So- 3 Pi- The pcr-axis kinetic energies [22j for this 
scan correspond to a horizontal (vertical) temperature 
Th(Tv) — 2 /iK (3 /iK), resulting in the t ypi cal thermal 
tail towards the red side of the resonance [27[. The solid 
line is a result of solving Q 



p(x,t) = -2K in (x 



p 2 (x,t) p{x,t) 



(3) 



with single particle density p, thermally-averaged in- 
elastic collision rate K m {x) = (K^k, A, x)) [22|, and 
vacuum- limited trap lifetime Tb g ~ 1.3 s. Equation [3] is 
solved at each x and we integrate the density to obtain 
the number of atoms at the end of the PA pulse. We 
then use Eq. [3] to fit to the experimental data [22| and 
extract £ op t and the position of the line center. Figurc[3jD 
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FIG. 3: (a) Typical PA loss feature for n=-2 in the low inten- 
sity regime at 7 av = 7 mW/cm 2 , with density- profile-averaged 
PA intensity 7 av [53] . (b) Time evolution of the trapped 
sample with (circles) and without PA light (squares). The 
two-body loss curve with PA is fit with a thermally averaged 
model (solid curve), (c) Linear increase of l op % with 7 av for 
iopt <S l/(2(fc)). (d) Molecular line center shift Av for large 
7 a v and decreased tpa- For each n, OFR parameters from the 
coupled-channels calculation and the experiment are summa- 
rized in the table at the bottom. Here, v n is the zero- intensity 
molecular line center with respect to 1 So- 3 Pi, and Av/I av 
characterizes the molecular ac Stark shift 1231. 



shows the time dependence of the same photoassociative 
loss process. Two-body loss is evident under resonant PA 
light. 

From the experimental data we extract two indepen- 
dent quantities: l pt1m an d an increased molecular loss 
rate 7 ~ 2.7j m . We have ruled out magnetic field or 
PA laser noise as a source of broadening. Instead, we 
conclude that this extra broadening is related to a faster 
molecular deca y ra te, which is consistent with our earlier 
measurements [2?5j ] and Rb results [1] . 

The measurements were performed for a range of £ op t 
by adjusting / av . Multiple molecular resonances were 
measured and results for n=-2 are shown in Fig.[3h. The 
optical length data is fit with a linear function and the 
results are summarized in the table at the bottom. The 
fit coefficient lopt/I&v is given by the free-bound Franck- 
Condon factor and decreases drastically with decreasing 
n. Figure [3ji exemplifies similar measurements done to 
determine the line shift Av with 7 av . Linear shift coef- 
ficients A^// av and zero intensity line positions v n with 
respect to the atomic transition are also shown in the 
table. The sign and magnitude of Ai//I av are consistent 
with the predictions in Rcf. 23]. 

At larger optical lengths (£ Q pt7m/7 ~ 100an)j clastic 
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collisions start to influence the dynamics of the system. 
We show the atom loss with respect to PA laser detuning 
for n=-2 in Fig. Both in-trap size and kinetic energy 
are measured by absorption imaging [22j]. Far detuned 
from the resonance, the gas is almost ideal, as shown 
by the persistent kinetic energy inhomogeneity along H 
and V in Fig. [4Jd. On resonance (vertical dashed line), 
inelastic collisions dominate and cause heating. For red 
detuning from the molecular resonance, the temperatures 
approach each other by cross-thermalization [29[. 

The measured cloud widths wh and wy confirm that 
the potential energy follows the kinetic energy (Fig. Hk) 
since particles oscillate in the trap many times between 
collisions. Similar measurements were performed for 
n=-3 and n=-4, and we find that the same dispersive 
behavior in temperatures and widths appears around 
2(fc)-^opt7m/7 ~ 30% at tpa = 200 ms. The data can 
be understood by a simple picture of competition be- 
tween K e \ and K[ n , which average differently with k in 
a thermal sample [see Eq. (JTJ], and thus peak at differ- 
ent values of A. Elastic collisions cause cross-dimensional 
thermalization and tend to equalize Th and Ty . Since in- 
elastic collisions predominantly remove cold atoms from 
the densest part of the cloud, the resulting loss increases 
the average system energy via anti-evaporation. 

This behavior is confirmed by a Monte-Carlo simula- 
tion, where 55 x 10 3 particles are simulated and each 
particle undergoes elastic and inelastic collisions with an 
initial phase-space distribution matched to the experi- 
mental conditions 



22] 



The solid lines overlaid on the 
experimental data in Fig.|4]are the simulation results. An 
average ratio of elastic to inelastic collisions per particle 
from the simulation is shown in Fig. HJi. The disper- 
sive shapes are also predicted by the coupled-channels 
model (see Fig. [2b) and their shape is sensitive to 7. 
Combined with the low £ op t data in Fig. [3l the entire 
simulation reproduces the experimental data only when 
7 = 2ir x 40(5) kHz without other free parameters. 

We conclude that the isolated resonance approxima- 
tion universally describes OFR in the vicinity of each 
resonance. The coupled-channels calculation includes all 
interference effects between resonances, and differs from 
the isolated resonance approximation at large detuning 
between resonances. Our experiment contradicts previ- 
ous predictions based on extrapolations of an isolated 
resonance to large detunings [id l20j . We have validated 
the linear line strength scaling and linear resonance shift 
with I and have observed a clear modification of both 
a(k) and b(k). For the values of i ptlm/j achieved here, 
inelastic losses still contribute significantly and (K e \/Ki n ) 
becomes even less favorable with decreasing T. However, 
the OFR effect can modify interactions in a degenerate 
gas of alkaline earth atoms and the desired change of 
a(k) is achieved at the smallest A/7 constrained by both 
molecular and atomic loss processes over a given experi- 
mental timescale 12211 . 
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FIG. 4: Elastic contribution to the scattering cross section 
for n=-2 at J a v = 22 mW/cm 2 (open circles) and results of a 
Monte-Carlo simulation (solid lines) using Eq. [T]in a crossed 
dipole trap for £ pt7m/7 = 140an. (a) Atom loss as a function 
of PA laser detuning from the atomic 1 So- 3 Pi resonance. In 
panels b and c, blue (red) data points and solid lines indi- 
cate the corresponding quantities for the vertical (horizontal) 
trap axis, (b) Change in kinetic energy derived from time-of- 
flight images, and (c) potential energy change corresponding 
to varying in-trap density profile, (d) The resultant ratio of 
elastic and inelastic collisions per particle, averaged over tpa- 
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SUPPLEMENTARY MATERIAL 
Scattering matrix, cross sections, and collision rates 

The collision of two identical bosons in the low energy 
near-threshold limit can be described in a scattering ma- 
trix treatment by a single s-wave scattering matrix ele- 
ment, S(k) = e 2lTl ^ k \ where hk is the relative momen- 
tum of the collision pair and rj(k) is the scattering phase 
shift due to interactions. When there is loss of scattering 
flux from the entrance channel, as is the case of a decay- 
ing optical Feshbach resonance, it is convenient to define 
a complex energy-dependent scattering length a(k), re- 
lated to the complex phase rj(k) and S-matrix element 
by 0,1 



a{k) = a(k) — ib(k) 



tanr?(fc) _ 1 1 - S(k) 
k ~~ ikl + S(k)' 



(4) 



This reduces to the usual definition of complex scattering 
length as k — > 0. The respective elastic and inelastic s- 
wave collision cross sections are 0, @| 



<T cl (k) = ^\l-S(k)\ 2 = 8n\a(k)\ 2 f(k), 



a in (k) 



2 -l(l-\S{k)\ 2 ) = ^b{k)f{k), 



(5) 



where f(k) = [1 + fc 2 |a(fc)| 2 + 2kb(k)]~ 1 -> 1 as k -> 
and \i is the reduced mass of the pair. The second 
relations in Eq. [5] follow from the definition in Eq. |4] 
The corresponding collision rate coefficients are 



Ad(fc) = — cr e i(fc), 

hk 

K ilL (k) = — tr in (fc). 



(6) 



When the entrance channel is coupled by a single fre- 
quency laser to a single excited molecular bound state 
of the pair, the isolated resonance approximation for the 
the field-dressed s-wave scattering matrix element can be 
written as 0, [H[ , 



S(jfc) ~ e 2lVb ^ k) [l-i 



r s (fc) 



E/H+ A + i(7 + r s (fc))/2 

' (7) 

Here E = h 2 k 2 is the collision energy with fi = ms r /2. 
The laser detuning A from the PA resonance includes the 
light shift Q to simplify the presentation (see Eq. [201 for 
a full expression). 

The stimulated width T s (k) = 2k£ opt (k)j m is the in- 
duced coupling between the free particle state \E) and 
the excited molecular state Here 7 m = 2~f a is the 
linewidth of the molecular transition and 7 a = 2tt x 
7.5 kHz is the atomic linewidth. We have allowed for 
extra molecular losses by letting 7 > 7 m . The optical 
length 



'-opt 



\{n\E)\ 2 ^ 
16nc k 



(8) 



where c is the speed of light, is a PA resonance- dependent 
line strength parameter [5j, [Z( • Under the Wigner thresh- 
old law, the free-bound Franck-Condon factor per unit 
collision energy |(n|_B)| 2 oc k, and we checked numeri- 
cally that it is a very good approximation to take £ opt 
to be independent of collision energy for temperatures 
< 10 /iK. The optical length is also independent of the 
atomic oscillator strength, and it scales linearly with PA 
intensity /. 

Combining Eq.[7]with Eqs.|4]|6]gives the near-threshold 
isolated resonance approximation to the OFR inelastic 
loss rate coefficient: 



K in {k) 



4nh 



M (A + E/h) 2 / 1 2 +[l + 2k^f^} 2 /4 



(9) 



The background s-wave scattering length a Dg is defined 
in the absence of a resonance as 



lim tan? ? b s( fc ) 



k 



(10) 



If we neglect a Dg for 88 Sr, we find a simple expression for 
the OFR-induccd clastic collision rate coefficient for an 
ideal ultracold gas 



K el (k) ~ 2k 



•'opt'Tr? 



K ia (k). 
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Power broadening is included via T s (k) in the denom- 



inator of Eq. which corresponds to 2k 



the 



denominator of Eq. |U From the expressions for elas- 
tic and inelastic collision rates, we see that the relevant 
strength parameter in the presence of extra molecular 
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loss is the rescaled optical length ^ pt-y-- Elastic colli- 
sions and power broadening only become important for 
the dynamics when 

2fc 4pt7m _ L (12) 



Coupled Channels Calculations 

The simplest way to do a calculation of the elastic and 
inelastic collision properties of an optical Feshbach res- 
onance is to set up a coupled channels model for the 
collision of field-dressed states in a laser field of fixed fre- 
quency 0, 11 • We use the simplest three-channel model 
that is sufficiently complete to represent nonpcrturba- 
tively the optically induced S(k) for an interfering spec- 
trum of excited state resonances. One channel represents 
the ground state with potential V g (R) and two channels 
represent the two excited states 0„ and l u that corre- 
late with the separated atom limit 1 So + 3 Pi with re- 
spective complex potentials Vq u and V\ u . These poten- 
tials are shifted to include the asymptotic detuning of 
the laser frequency from atomic resonance. It is neces- 
sary to include both U and l u excited states to correctly 
calculate the complex scattering length near an isolated 
resonance of either state. This is because the molecular 
J = I excited states (J= total angular momentum quan- 
tum number) must become a mixture of s and <i-waves 
at long range to properly go to separated atoms with 3 Pi 
quantized in a space frame instead of a molecular body 
frame Q ; otherwise, spurious 1/i? 3 resonant dipole terms 
would get mixed into the ground state potential and give 
invalid threshold scattering lengths. 

The coupled channel potential matrix V(R) is found 
using the asymptotic representation of the two excited 
J = I molecular channels in terms of a pair of excited 
state s- and <i-waves for the e parity block in Table 1 of 
Reference [(§: 

/ Vg Vo P t \ 

V= \ V opt \{V 0u +2V lu ) ^(V lu -V 0u ) , 
V &(V lu -Vou) l(2V 0u + V lu +QB) J 

(13) 

where 6B = 6ft 2 /(2/ii? 2 ) is the <i-wave centrifugal en- 
ergy. Here the row labels, in order from top to bottom, 
represent the \j£JM) = |0000) ground state and 1 1010) 
and 1 1210) excited state channels of Ref. where j is 
the atomic angular momentum for the collision partner 
to the 1 So atom, i is partial wave, and JM represent 
the total angular momentum and its projection for the 
pair of interacting atoms. Since the excited bound states 
are relatively short range, we use the nonretarded opti- 
cal coupling V op t = (27ri / 'c) 1 / 2 ^™, where d rn = \[2d,A 
and the atomic transition dipole oIa — 0.086816 atomic 
units corresponds to an atomic lifetime of 21.46 /is. A 



complete theory would use the retarded optical coupling. 
The zero of energy for ground state collisions is taken to 
be the lowest eigenvalue of this matrix to account for the 
energy shift of the field dressed molecule. The numerical 
S'-matrix is calculated using standard coupled channels 
methods to calculate the wave function for the coupling 
matrix in Eq. 1131 

While Reference Q simulates the spontaneous radia- 
tive decay of the excited state by introducing artificial 
channels, here we include a complex term —ihq/2 in each 
of these excited state potentials to simulate this decay, 
where in our numerical calculations we take 7 = 7,„. 
This leads to a non-unitary ground state S'-matrix ele- 
ment such that < 1 - \S{k)\ 2 < 1 gives probability 
of atom loss during a collision. In order to avoid prob- 
lems with spurious asymptotic decay of the field dressed 
states, a tapered cutoff function is introduced so that 
decay is only turned on inside some characteristic dis- 
tance, which we take to be 500 ao- We have verified that 
a(k) is independent of this choice over the full range of 
detunings, even between resonances, and b{k) is indepen- 
dent of this choice out to many linewidths away from an 
isolated resonance. However, in the far wings of a reso- 
nance between isolated resonances, b{k) will depend on 
this choice, although molecular loss tends to be small in 
these regions. This dependence on cutoff is because our 
non-unitary theory does not distinguish between atomic 
and molecular light scattering when the atoms are far 
apart; the cutoff selects which distance inside of which 
an excitation is considered to be a molecular loss process 
instead of an atomic one. Thus the cutoff is taken to be 
on the order of but larger than the outer turning point 
of the excited molecular bound states. 

Finally, we ignore coupling of the excited J = 1 molec- 
ular levels back to d-waves in the ground state. This 
means that our model only gets the s-wave but not the 
d-wave contributions to the excited state light shifts [||. 
It would be straightforward to add an extra ground state 
channel to the numerical calculation to account for this. 

While there are a number of ways our simple model 
can be improved, it gives a practical way to calculate 
the needed complex scattering length for characterizing 
the intercombination line OFR spectrum of alkaline earth 
atoms. The model gives two key results. First, it should 
correctly describe the variation of the real part of the 
scattering length as laser frequency is tuned across mul- 
tiple resonances. Second, the isolated resonance approx- 
imation should correctly describe the molecular loss rate 
out to order of hundreds of line widths away from indi- 
vidual isolated resonances. 



Experimental Setup 

Details of the experimental setup can be found in 
Ref. ||. Atoms are loaded into an optical dipole trap 
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from a magneto-optical trap operating on the 1 So- 3 Pi 
intercombination transition. A horizontal (H) and a ver- 
tical (V) beam intercept in the x — z plane to form a 
crossed dipole trap at the origin shown in Fig. [5] (draw- 
ing is to scale). Both beams are derived from the same 
1064 nm laser and are linearly polarized along the y axis 
with 1/e 2 waists 63 /tm and 53 /mi, respectively. The 
optical beams point in the directions 

kn = — cos8hx — sin 6*#i, 
ky = — cos 6yz + sin 8yX, (14) 
fcpA = — cos^paA + sin^pAy, 

with Oh = 16.0°, 8 V = 14.4°, and 6 PA = tt/8. The beams 
are linearly polarized along e# = ey = y and epA = z. 
The image plane is spanned by cos ^ x — sin ^y and z. 
A bias magnetic field of B z ~ 100 mGauss defines the 
atomic quantization axis. 




FIG. 5: Geometry of the experiment in the absorption image 
plane. The directions x, y, and z define the lab frame, where 
both gravity g and bias magnetic field B are parallel to z. 
Symbols k are beam directions, and e are beam polarization 
vectors, where subscripts H , V , and PA indicate horizontal, 
vertical, and PA beams. H and V Gaussian beam profiles 
are shown in blue outline, PA Gaussian beam profile in red 
outline. 

The focal condition is not critical to model the poten- 
tial sufficiently well because of the large Rayleigh ranges 
(> 8 mm), and the Gaussian beam isosurfaces at the 
beam waist are almost cylindrical (blue outlines). Thus, 
the dipole trap potential due to each beam is modeled as 

Ui(x) = -t/^exp j--^ ([e.i ■ x} 2 + [(ki x ii) ■ ar] 2 ) j , 

(15) 

with individual Gaussian beam trap depth [9[ 
p. 

U % T = ^jRe ai Sn (1064 nm), (16) 

and total laser power Pi. The real part of the dynamic 
1 So polarizability at 1064 nm is Re ais o (1064 nm) ~ 
239 a.u. fioj ]. with atomic unit of polarizability 47reoag. 



The full model potential includes the gravitational ac- 
celeration g pointing along — z and is given by 

U(x) = U H (x) + U v (x) + m Sx gz. (17) 

Gravity sets the trap depths of ~7 /iK and ~T5 /iK along 
V and H. The two graphs on the right hand side of Fig. [5] 
show cuts through the model potential, which has been 
adjusted to match the measured trapping frequencies. 
An isosurface (dark blue) of the model potential at 7 /iK 
is shown in the zoomed-out portion. 

Typical kinetic energies are 2-4 /tK, and in-trap Full 
Widths at Half Maximum (FWHM) are 45-55 /tm. The 
PA beam (red outline) propagates along the horizontal 
image axis in the x — y plane with a waist iupa = 41 /tm 
and is linearly polarized along z. Although the PA beam 
diameter 2v21n2 x wy>a — 97 (im is larger than the typ- 
ical cloud FWHM, we use a density-averaged intensity 
7 av = J d 3 xp(x)I(x)/ J d 3 xp(x) to characterize the PA 
intensity interacting with the atoms. Typical values are 
I av ~ (0.6 — 0.7) x /pk, where I p k = is the Gaussian 

beam peak intensity for total beam power P. 

In addition to scattering from the atomic transition 
with rate r sc , the PA beam also adds to the optical trap 
slightly. For large intensities and small u>pa and espe- 
cially in a standing wave configuration, this effect can be- 
come important. In this work, the additional trap depth 
introduced by the PA beam is typically < 0.1 /iK. 

Imaging Procedure and Analysis 

We use a two-component long- working-distance micro- 
scope with 3 /im resolution at 532 nm. The imaging lens 
is placed outside the vacuum chamber at a distance of 
150 mm from the atomic cloud. The second part of the 
telescope is an eyepiece mounted to a CCD camera with 
1024 x 1024 square pixels of 13 /tm width. The eyepiece 
includes an interference filter at 461 nm (intensity trans- 
mission 0.45 within 10 nm bandwidth, ND 5.0 otherwise). 
Time of flight measurements were used to determine an 
imaging magnification of 6.5, corresponding to 2.0 /tm 
per pixel. We tested the imaging system resolution in 
a separate test setup that included the effect of an anti- 
reflection (AR) coated vacuum viewport, same as the one 
used on our vacuum chamber. Using incoherent white 
light filtered by the blue interference filter with transmis- 
sive test targets, resolutions of ~ 4 /tm were measured. 
Absorption imaging with coherent light causes etaloning 
between the CCD chip and the glass plate that covers it. 
To reduce fringing, the glass plate surface is wedged at 
~ 1° and AR coated for the imaging wavelength. 

Without the elastic contribution of the OFR effect, 
the 88 Sr sample does not thermalize on experimental 
timcscales due to its small <2b g , and it is important to 
obtain information about both the potential and kinetic 
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energies of the gas. We interleave experimental runs be- 
tween in-situ imaging of the sample and imaging after the 
atoms were dropped for 1.5 ms by switching off the opti- 
cal dipole trap. The first image gives information about 
the spatial density of the atomic cloud, and the second 
image measures the kinetic energy by time-of-flight ex- 
pansion. We typically repeat this interleaved sequence 
3-5 times and average the resulting in-situ and time-of- 
flight absorption images before performing data analysis. 
We also calculate the pixel-by-pixel studentized standard 
deviations which are used as weights in fitting the mean 
image and also enter the fit parameter error estimates. 
Both the mean and standard-deviation images are then 
rotated (with bilinear interpolation) by 14° about the 
center of the atom cloud to transform the images into 
the eigenframe of the trap. The pictures are then inte- 
grated along the horizontal (vertical) axis and are fit to 
extract wh (w) with one-dimensional Gaussian distri- 
butions including a background offset and linear slope to 
account for residual fringing and CCD readout noise in 
the absorption image. 

A common problem in imaging atomic clouds with 
large optical depth is light that does not interact with 
atoms hitting the camera in the same position. Common 
causes arc forward scattering of photons into the imag- 
ing path or mixed probe polarization. Probing the 
1 Pi transition is insensitive to probe polarization even 
under a small bias magnetic field. To identify the pos- 
sibility of a limiting optical depth OD co iii ng , a series of 
images were taken with increasing atomic density. We 
compared the peak optical depth in-situ versus the pixel- 
summed optical depth either of an in-situ or a corre- 
sponding time-of-flight image, and we estimate a con- 



servative limit OD c 



> 3.5. Typical in-situ optical 



depths in this paper were below 2, so that we can neglect 
the effect. 

For probe intensities larger than 0.1/ sa t, with satu- 
ration intensity I sat ~ 40 mW/cm 2 , saturation of the 
imaging transition becomes important. To account for 
the saturated absorption, we correct the full expression 
for the OD to the linear Beer-Lambert regime [1, [ill ]. 
For each picture the saturation correction is applied on 
a pixel-by-pixel basis. 

While the shadow image is forming and the atoms are 
scattering light, their momentum undergoes a random 
walk due to the spontaneous emission of blue photons. 
By examining the shadow image versus a reference im- 
age without atoms, we extract the number of photons 
that the atoms removed from the probe pulse. By vary- 
ing the probe time, we find the cloud expansion as a 
function of the number of photons scattered and apply a 
corresponding correction to the cloud widths from both 
in-situ images and TOF images [|| . 

We extract the integrated OD of the mean image to 
obtain the atom number. Corrections for probe heating 
are then applied to the fitted widths. The in-situ image 



widths are used to calculate the mean atomic density. Fi- 
nally, we calculate a time-of-flight temperature corrected 
for finite size effects by the in-situ widths. 



Loss Spectra (2fc£ op t7m/7 <S 1) 

To extract £ op t"fm from the loss spectra, data was 
fit to an approximate expression for the integral of the 
atomic density, where the density after the PA pulse was 
calculated via the differential equation p = —K ln p — 
p/ T bg [13] • The initial density distribution was modelled 
as Gaussian with wh and wy ■ The width into the plane 
of the images was extrapolated from a Monte-Carlo sim- 
ulation of the trap model potential. The fitting function 



N = 



du 



y/u e" 



1 + A^ in (A T) r cff p e 



(18) 

with normalized dimensionless trap length scale u. This 
expression describes the number signal normalized to the 
measured atom number after the PA pulse of duration 
tpa and when the detuning from molecular resonance 
is large. The quantity t c s = T bg e' rhold / Tb <5(e TpA / rb <5 - 1) 
accounts for the fact that atoms are lost to one-body 
decay during the PA pulse with lifetime Tb g and that we 
hold the atoms for Thoid = 100 ms between the end of 
the PA pulse and the imaging. The number po is the 
peak density after the PA pulse and at large detuning. 
The definition of all other symbols follows the main text. 
To make the fit numerically tractable, the integral was 
approximated as a ten-term sum using Gauss-Laguerre 
quadrature. 

The thermally-averaged inelastic collision rate per par- 
ticle K m was derived from the s-wave scattering matrix 
of Bohn and Julienne, using the definition of the optical 
length £ opt = T s /2kj m in terms of the stimulated width 
Ts [J,0]- Here, we use a thermally- averaged inelastic rate 
coefficient K m = (Hka- m /fj,), where a- m = 2p-(l — \S\ 2 ) 0]. 
The angular brackets denote a thermal average over the 
initial collision momenta. This average was performed 
on the inelastic rate coefficient (rather than the density 
after the PA pulse) under the assumption that prior to 
the PA pulse, many single-particle velocities exist in ev- 
ery differential volume in the trap. Since single collisions 
take place in a small volume, the spatial dependence of 
K- ln was carried through into the differential equation for 
the density rather than averaged out. 

To fit the inelastic loss spectra, K m was expressed as 



T) 



D 



8y / 7rMo P t7 n 



dr] 



D 2 



h 



r = 7 + 2k th e opt j my /rj, 



f r 2 /4' 

-U/UQ 

) 

(19) 
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with dimcnsionlcss relative momentum magnitude rj = 
k 2 /k^ h , thermal momentum Hkth = V^fJ-ksT, center-of- 
mass PA photon recoil energy hv Tec , trapping laser ac 
Stark shift at the center of the trap v s , and Planck's 
constant h = 2irh. The PA laser with optical frequency 
Vi is detuned from the PA resonance by 



A = 14 - [vfSo - 3 Pi) 



■Au(I)], 



(20) 



where f( So — 3 Pi) is the atomic transition frequency, 
and hv n is the energy of state n with respect to the free 
threshold. The detuning term accounts for the ac 

Stark shift of the molecular resonance with respect to /. 

The integral in Eq. [15] was approximated as a 53-term 
sum using Gauss-Laguerre quadrature. The quantities 
(A>pt7m)i T, uq, and a term added to the detuning to 
represent the line center were allowed to vary. The pa- 
rameter T is used as a check against the experimentally 
measured temperatures and agrees well with the experi- 
ment. The Stark shift term v s e~ u ^ u ° was included to ac- 
count for the broadening of the atomic loss profile to the 
blue side of a molecular resonance due to the ac Stark 
shift induced by the trap. The molecular line width 7 
was extracted from the Monte-Carlo simulation in the 
next Section. 



Monte-Carlo Simulation (2fc£ op t7?n/7 ~ 1) 

To model the thermodynamic effects caused by the in- 
terplay of elastic and inelastic collisions in an anharmonic 
trap with evaporation, we use a Monte-Carlo simula- 
tion since analytic expressions such as those presented 
in the previous Section are not available. The simula- 
tion is based on classical particles moving in a conser- 
vative model potential that includes the Gaussian beam 
shapes as well as gravity. A commonly used method [13- 
15| to include collisions in such simulations is due to 
Bird [lj§ and detailed discussions of the method in the 
context of ultracold atoms in optical traps can be found 
in Rcfs. 



171 Il8f . The simulation procedure here is de- 



scribed in more detail in Ref. |8[. 

The simulation uses the elastic and inelastic cross 
section formulas derived in Section "Scattering Matrix, 
Cross Sections, and Collision Rates" : 



<7in(fc) 



4tt 



cr c i(fc) ~ 2k 



k (A + E/h)y^ + (1 + 2*^)74' (21) 

n 



1 



'-cr in {k), 



where the second relation is an approximation valid in 
an ideal ultracold gas like 88 Sr where the background 
scattering length ab g can be neglected. For a gas at non- 
zero temperature, the detuning A includes the atomic 
motion and the trap ac Stark shift as in the previous 
section. 



The initial particle distribution is synthesized by drop- 
ping atoms into the model potential and letting them 
evolve for several hundred ms without collisions (using 
an embedded Rungc-Kutta method) . Each particle is ini- 
tially generated from independent Gaussian distributions 
along the trap eigenaxes, both in position and velocity. 
The initial widths of these Gaussian distributions are ad- 
justed until the particle distribution after having settled 
in the model potential matches the experimental in-situ 
and TOF data when the PA laser is far detuned from the 
OFR resonance. 



To calibrate the Monte-Carlo simulation, we modelled 
a three-dimensional isotropic harmonic trap of mK trap 
depth and checked the thermodynamic effects of inelastic 
and elastic collisions independently. Using a Gaussian 
initial phase-space density at several fj,K and only al- 
lowing inelastic collisions at a velocity-independent cross 
section, we recovered heating rates per particle that 
match the corresponding analytic expressions Simi- 
larly, we check cross-dimensional thermalization rates at 
a velocity-independent clastic cross section. We found 
that the average number of elastic collisions required 
for cross-dimensional thermalization in our simulation 
agrees [|| with the analytic expressions from Ref. [Hf. 

The Monte- Carlo simulation also includes the effect of 
scattering from the atomic line. For the largest optical 
lengths in Fig. 4 of the main paper, the atomic scattering 
rate is ~ 6 s^ 1 corresponding to 1.2 photons scattered per 
atom during the 200 ms exposure time. At these scatter- 
ing rates, the photon absorption along the PA laser di- 
rection and the random reemission change the resulting 
cloud widths by less than 10%. Due to the spontaneous 
reemission as a spherical wave, the mean free path of a 
scattered photon in a sample of OD ~ 2 is still larger 
than the typical cloud radii and we estimate that radia- 
tion trapping should only introduce a small correction. 



The final particle distribution is then imaged by ran- 
domly scattering photons off of each particle and forming 
a histogram from the positions of the scattering events 
over the image pulse time of 50 //s, either in the trap 
or after TOF expansion for 1.5 ms. Random rescatter- 
ing and corresponding velocity and position evolution are 
included, and the position histogram is blurred to ac- 
count for the 4 fim image resolution. The final image is 
then analyzed in the same way as the experimental data 
to extract in-trap widths and TOF temperatures. By 
adding the final measurement step instead of calculating 
the covariance matrices of the position and velocity dis- 
tributions directly, the quantitative agreement with the 
experiment was improved significantly in the regime of 
large inelastic losses. 
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Maximal scattering length modification 

In the following, we build on the understanding from 
our current experiment and try to estimate the maxi- 
mum scattering length under the most ideal conditions. 
We assume that atomic scattering is the dominant loss 
mechanism, that k — ► 0, and that there is no extra molec- 
ular loss, such that 7 = 2 r y a . There is also no spatial 
inhomogencity, and thus it docs not reflect the current 
experimental conditions. Theoretical calculations of the 
molecular line strength factors are used to make these 
estimates. Also note that the line strength factors are 
only order-of- magnitude estimates beyond n < — 5. 
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The Optically-Modified Scattering Length 
Constrained by Atomic Light Scatter 



FIG. 6: Equation [26] plotted near molecular resonance. Here 
an £ corresponding to the —23.264 MHz line has been used, 
T sc = 1 s"\ 7a = 2tt x 7.5 kHz, and 7 sat = 3 fiW/cm 2 . This 
behavior is well approximated by Aadis of Eq. 1271 



In the k — > limit, the scattering length change Aa 
due to an optical Fcshbach resonance is given by 



^-opt 



(S - (5 )7r, 



(<5-<5o) 2 +7^/4 



21 



opt 



(6 - S Q )j a 

(s - s y 



ll 



(22) 

where S is 2ir times the detuning from atomic resonance, 
Sq is the difference between the molecular and atomic 
resonance frequencies, j m is the molecular linewidth, j a 
is the atomic linewidth, and y m = 2^/ a . The quan- 
tity £ opt = £1, where / is the PA laser intensity and 
\( n \E)\ j s a rcsona nce-specific constant that de- 
scribes the Sr+Sr molecular structure. 

When using an OFR, 8 must be chosen such that in- 
elastic collisional loss is small. In this regime, loss due to 
atomic light scattering dominates. The atomic scattering 
rate T sc is given by 



r — 

J- sc — 



la 



So 



2 1 + s + 4(<V 7a ) 5 



(23) 



Here sq = I/I sa ,t, where J sa t is the saturation intensity 
of the atomic transition. Typically the experiment de- 
termines the maximum allowable value of r sc , thereby 
constraining /. Under this constraint, / can be expressed 
as 



I = I» 



2r sc /7 a 

1 - 2F sc / 7a 



7a 2 



8/ sa 



7a 



(24) 

where it has been assumed that j a ^> T sc and 5 ^> j a , the 
latter of which is true due to strontium's narrow inter- 
combination line that is used in our work. The maximum 
£ opt for a given detuning is simply 



V = = 8£/ sat — — 



7a 



la 



(25) 
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Detuning (GHz) 

FIG. 7: Equation [26] plotted for large detunings. This plot 
was generated with the same values for £, <5o, F sc , y a , and 
/sat as those used in Fig. [6] This behavior of Aa for large 
detunings is well approximated by Aan n of Eq. 1281 As will 
be discussed in Section "The Large Detuning Case," these 
detunings are unphysical given the molecular structure and 
the coupled-channnels theory discussed in the main text. 



Thus Aa, constrained to a certain r sc , is 



Aa = 16£J 8at — 

la 



S 2 (5 - S h a 

l 2 a(S-5 ) 2 +l? L 



(26) 



This expression is plotted in Figs. [5] and [7] Within tens 
of linewidths 7 m of 6 = Sq, the factor 5 2 in the numerator 
of Eq. [26] is slowly varying and can be approximated as 
5q. Let Aadis refer to Aa in this regime. Therefore, 



Aa dis = 16£I sat — 



Si (6-<5o)7a 

la J ll (<5-5 ) 2 + 7 2 : 



(27) 



which has the dispersion shape of Fig. [5] When <5 ^> <5q, 
|Ao| varies linearly with 5. Let Aaij n be the expression 
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for Aa when S 3> 5q . 

Aaii n - 16£/ sat f £^ , (28) 

which describes the linear behavior evident in Fig. [7] 

The Maximum Useful Scattering Length Near a 
Molecular Resonance 




-23.34 -23.3 -23.26 -23.22 -23.18 
Detuning (MHz) 

FIG. 8: Equation [26] plotted in blue and Eq. [29] plotted in 
red. This plot was generated with the same values for f, <5o, 
r sc , 7a, and J sa t as those used in Fig. [6] Note that K- ln is 
significant near a molecular resonance but that this quantity 
drops off faster than Aa. 

Near a molecular resonance where |Aa| has dispersive 
behavior, it follows from Eq. [27] that |Aa| has a local 
maximum at \S — Sq\ ~ ±7o- However, at these small 
detunings inelastic collisional loss must be considered. 
For k — > 0, the inelastic rate coefficient K- m is given by 

= 47rM opt 7rn 

A 4 (S-6 r +1 l/4 
_ 16TTM op t ~fl 

/x {S- So) 2 + i 2 a ■ 

Putting Eq. [25] into Eq. [29] it can easily be shown that 
the maximum value of K m , denoted by K™ ax , occurs at 
S rj So, consistent with the case when £ opt in Eq. |2"21 is 
not constrained by atomic light scatter. 

As Fig. [8] shows, inelastic loss (described by K- m ) is 
significant near So- Nevertheless, this figure makes it ap- 
parent that one can choose a detuning that is just outside 
of the influence of K- m but that still yields a significant 
Aa. This detuning is implied upon setting K m = kA™, 
where k describes the fraction of molecular loss deemed 
acceptable by the experiment. Solving this equation for 
detuning yields that S—5 ~ ±7a/ Evaluating Eq.l26l 
at this detuning produces Aa ma x, 

Aa max « ±16£T sa tV^ (^j , (30) 



n 


6 (GHz) 


<5lin (GHz) 


-2 


-0.023 


-7.22 


-3 


-0.221 


-655 


-4 


-1.084 


-1.57 xlO 4 


-5 


-3.460 


-1.60 xlO 5 


-6 


-8.400 


-9.40 xlO 5 


-7 


-17.0 


-3.84 xlO 6 



TABLE I: Theoretical values for So for the n = —2 through 
-7 U resonances and the associated values for 5u n - Here n = 
0.01. 



where Aa max is the maximum scattering length change 
when Aa is dispersive and inelastic loss and atomic light 
scatter are negligible. The intensity corresponding to 
Aa max , given by Eq. is 

/Aa -" 8/sat (^)0 2 - (3i) 



The Large-Detuning Case 

Although |Aa| in Eq. [26] seems to increase without 
bound when S ^> So, coupled-channels OFR theory 
(which so far has not been considered in this discussion) 
dictates that if |Aa| is based on the nth resonance, the 
modification to the scattering length will vanish at some 
point before S is equal to the detuning of the (n — l)th 
resonance. See the main text for a further discussion of 
this coupled- channels effect. One must consider whether 
the molecular detuning 5 — So can be sufficiently large to 
take advantage of the linear increase with S but not so 
far that coupled-channels effects become a concern. 

Let S\i n be the detuning at which the magnitude of 
Aaii n is equal to |Aa max |. It follows that |<5i; n | = 
V / K ( 5o/7a - ^ ^im is comfortably outside the regime where 
coupled-channels effects are a concern, then the maxi- 
mum useful scattering length occurs when one is many 
lincwidths 7 m detuned from So but not far enough de- 
tuned to be close to other resonances. However, theoret- 
ical values for So make it clear that detuning to Sn n would 
always require crossing another molecular resonance, so 
coupled-channels effects dictate that the maximum scat- 
tering length one can achieve is given by Aa max (Eq. [30"]) . 
occurring when S is as close as possible to a molecular 
line, while molecular losses are constrained to a given 
level. Table U considers 6 different molecular resonances 
and their associated values of Sn n , which are many orders 
of magnitude larger than <5q for subsequent resonances. 
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n 


5 (GHz) 


€ (a /(W/cm 2 )) 


Aa max 


(ao) (W/cm 2 ) 


-2 


-0.023 


6110 


5.97 


4.89 xlO -3 


-3 


-0.221 


32.8 


2.92 


0.443 


-4 


-1.084 


27.9 


59.3 


10.6 


-5 


-3.460 


3.30 


71.5 


108 


-6 


-8.400 


0.272 


34.7 


637 


-7 


-17.0 


0.012 


6.21 


2.61 xlO 3 



More significant gains can be achieved by shortening ex- 
perimental time scales to allow for larger values of r sc 
since |Aa max | is linear in r sc . 



TABLE II: Theoretical values for So and £ for the n — — 2 
through -7 U resonances are used to calculate the maximum 
useful scattering length. The intensities required for these 
scattering lengths are also included. The optimal detuning 
for each resonance is given by \8 — 5o\ — j a / 1 = IO7 for the 
value k — 0.01 that was used in the table. 



Conclusion 

In light of Section "The Large-Detuning Case" , we con- 
clude that for a given molecular line, the best optically- 
modified scattering length occurs for a laser detuned 
just outside the influence of inelastic loss and with an 
intensity just low enough to prevent heating due to 
atomic light scatter. This optimal detuning is given by 
S — 5 = 7 Q l-JH. Table ILT1 lists theoretical 0„ values for 
£ and Sq as well the associated values for Aa max and 
JAn m „. Values of r sc = 1 s _1 and k = 0.01 have been 
used. For the strontium 1 So- 3 Pi line, j a = 2ir x 7.5 kHz 
and J S at = 3 fiW/cm 2 . These numbers yield an opti- 
mal detuning of S — 5q = 27r x 75 kHz. The intensity 
column was included to point out that many resonances 
cannot be used on technical grounds (for instance in- 
tensities greater than a few kW/cm 2 are impractical for 
external cavity diode lasers). 

Under the ideal conditions of zero temperature, no spa- 
tial inhomogeneity, no extra molecular loss and the as- 
sumption that atomic light scattering is the dominant 
loss mechanism, Table |H] shows that the resonance loca- 
tion of the n = —5 line is the best compromise between a 
large detuning from atomic resonance (which allows for 
more PA laser intensity without atomic light scatter) and 
the fact that £ rapidly decreases with increasing Ac- 
cording to Eq. [3D1 only moderate gains in Aa max can be 
achieved by allowing for more inelastic loss (when detun- 
ing closer to a molecular line) since any increase in K ul 
will only be accompanied by a \/K ln increase in | Aa max | . 
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